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ABSTRACT 

We present a simple spin evolution model that predicts that rapidly rotating accreting neutron stars 
will mainly be confined to a narrow range of spin- frequencies; P — 1.5—5 ms. This is in agreement 
with current observations of both neutron stars in the Low-Mass X-ray Binaries and millisecond radio 
pulsars. The main ingredients in the model are: i) the instability of r-modes above a critical spin rate, 
ii) thermal runaway due to heat released as viscous damping mechanisms counteract the r-mode growth, 
and iii) a revised estimate of the strength of dissipation due to the presence of a viscous boundary layer 
at the base of the crust in an old and relatively cold neutron star. We discuss the gravitational waves 
that are radiated during the brief r-mode driven spin-down phase. We also briefly touch on how the new 
estimates affects the predicted initial spin periods of hot young neutron stars. 

Subject headings: dense matter - gravitation - stars: neutron - stars: rotation — stars: oscillations 



1. INTRODUCTION 



Kokkotas (2000 ) for a review) to dissipate the accreted 



The launch of the Rossi X-ray Timing Explorer (RXTE) 
in 1995 heralded a new era in our understanding of neu- 
tron star physics. Detailed observations of quasiperiodic 
phenomena at kHz frequencies in more than a dozen Low- 



angular momentum from the neutron star 

In this Letter we reexamine the idea that gravitational 
waves from unstable r-modes provide the agent that bal- 
ances the accretion torque. This possibility was first an 



Mass X-ray Binaries (LMXB) strongly suggest that these 



jyatomQ contain rapidly opinning ne utron otaro (for a re 



alyzed in detail by Andersson, Kokkotas fc Stergioulas 
(199S| ) (but see also jBildstcn (1998| )). Originally, it was 



cent review, see yan der Klis (2000 )), providing support 
for the standard model for the formation of millisecond 



thought that an accreting star in which the r-modes 
were excited to a significant level would reach a s pin- 



pulsars (MSP) via spin-up due to accretion. 



to be answered by further observations and/or theoretical 
modeling. For example, we still do not know the reason 
for the apparent lack of radio pulsars at shorter periods 
than the 1.56 ms of PSR1937+21 (for a review of recent 
progress in the m odelling of rotating neutron stars, see 
Stergioulas (1998] )). The recent RXTE observations pro- 



equilibrium, very much in the vein of suggestio ns by Pa 

n"" -T"^ ' "j " H i " "£" iv r" i ""' . paloizou & Pringle (1971) and IWagoner (1984). Should 
Doupito tnoao advanccf] Dovoral dimciilt quoutionu remain , . — : : — a ■ 1 J — tttttt; ' r 3 , 

this happen, the neutron stars m LMXB would be prime 

sources for detectab le gravitatio nal w aves. Howeve r, as 
was pointed out by Levin (1999 ) and Spruit (1998 ), the 
original idea is not viable since, in addition to generat- 
ing gravitational waves that dissipate angular momentum 
from the system, the r-modes will heat the star up (via 
the shear viscosity that counteracts the r-mode at the rel- 
evant temperatures). Since the shear viscosity gets weaker 
as the temperature increases, the mode-heating triggers a 
thermal runaway and in a few months the r-mode would 
spin an accreting neutron star down to a rather low rota- 
tion rate. Essentially, this conclusion rules out the r-modes 
in galactic LMXB as a source of detectable gravitational 
waves, since they will only radiate for a tiny fraction of 
the systems lifetime. 

Other recent results would (at first sight) seem to em- 
phasize the co nclusion that the r-modes are not relevant 
for the LMXB. |Bildsten fc Ushomirsky (2000| ) investigated 
the effect that the presence of a solid crust would have on 
the r-mode oscillations. They estimated that the dissipa- 
tion associated with a viscous boundary layer that arises 
at the base of the solid crust in a relatively cold neutron 



vide a further challenge for theorists. Various models sug- 
gest that the neutron stars in L MXB spin rapidly, perhaps 
in the narrow range 260-590 Hz ( |van der Klis, 200C| ). Three 
different models have been proposed to explain this sur- 
prisin g result. The first model (due to White fc Zhaiip] 



(1997 )) is based on the standard magnetosphere model for 



accretion induced spin-up, while the remaining two mod- 
els are rather different, both being based on the idea that 
gravitational radiation balances the accretion torque. In 



the first su ch model for the LMXB (proposed by Bild ■ 
sten (199c|) and recently refined by Ushomirsky, Cutlci 



fc Bildstcn (2000| )), the gravitational waves are due to a 



quadrupole deformation induced in the deep neutron star 
crust because of accretion generated temperature gradi- 
ents. The second gravitational- wave model relies on the 
recently discovered r-mode instability (see Andersson & 
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star would greatly exceed that of the standard shear vis- 
cosity. Thus, Bildsten and Ushomirsky concluded that the 
r-mode instability would only be relevant for very high ro- 
tation rates, and could therefore not play a role in the 
LMXB. 

We have reassessed the effect of the viscous bou ndary 



layer (correcting an erring factor in the estimates of Bild 



sten & : Ushomirsky (200C)). Our new estimates show that 
the presence of the crust is important, but that the insta- 
bility operates at significantly lower spin rates than sug- 
gested by Bildsten and Ushomirsky. Once we combine our 
estimates with the thermal runaway (now due to heating 
caused mainly by the presence of the viscous boundary 
layer), that results as the star is spun up to the point 
at which the instability sets in, we arrive at a model for 
the spin-evolution of accreting neutron stars. Remarkably, 
this simple model agrees well with existing observations of 
rapidly rotating neutron stars, covering both the LMXB 
and MSP populations. 

2. DISSIPATION DUE TO A VISCOUS BOUNDARY LAYER 

The r-mode instability follows after a tug of war be- 
tween (mainly current multipole) gravitational radiation 
that drives the mode and various dissipation mechanisms 
that counteract the fluid motion. In the simplest model, 
the mode is dominated by shear viscosity at low tempera- 
tures while bulk viscosity may suppress the mode at high 
temperatures. At intermediate temperatures, the r-mode 
sets an upper limit on th e neutron star spin rate. In an 
interesting recent paper, Bildsten fc Ushomirsky (2000 ) 



estimate the strength of dissipation due to the solid crust 
of an old neutron star, and find that the presence of a 
boundary layer at the base of the crust leads to a very 
strong damping of the r-modes. 

While we agree with the main idea and the various as- 
sump tions made bv Bildsten and Ushomirsky. we would 



like t< p point out one important difference between their 
results and ones used previously in the literature. Their 
assumed timescale for gravitational radiation reaction dif- 
fers significantly from, for example, the unifor m density 
result derived by Kokkotas fc Stergioulas (1999 ) (and sub 



sequently used by several authors, see Andersson & Kokko- 
tas (200C| )). This is surprising since the uniform density 



result, which can be written 
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(where the negative sign indicates that the mode is unsta- 
ble) has been shown to be close (within a factor of two) 
to the results for n = 1 polytropes. M, R, and P repre- 
sent the mass, radi us and spin period of the star , respec- 
tively. In contrast, Bildsten fc Ushomirsky (200C| ) use the 
n = 1 polytrope result and argue that it corresponds to 
tg W ~ —146 s for a canonical neutron star rotating with 
a period of 1 ms, i.e. assume that radiation reaction is 
almost one order of magnitude weaker than in (|l|). This 
difference occurs because Bildsten and Ushomirsky have 
only rescaled the fiducial rotation frequency f^o = \J~xGp 
(where p represents the avera ge density) in terms of which 
the n = 1 polytrope results of Owen et al. (1998) were ex- 
pressed (t gw ps — 3.26(Qo/£l) e for a specific polytropic stel- 
lar model). Unfortunately, this procedure is not correct. 



From the fundamental relations, e.g. the formula for the 
gravitational-wave energy radiated via the current multi- 
poles, one can see that the gravitational-wave timescale 
should scale with M, R and P in the way manifested in 
(0). Thus, we believe that Bildsten and Ushomirsky un- 
derestimate the strength of radiation reaction significantly, 
which motivates us to reassess the relevance of the viscous 
boundary layer. 

We should, of course, emphasize at this point that our 
current understanding of the r-mode instability is based 
on crude estimates of the various timescales. In order to 
understand the role of the instability in an astrophysical 
context we must improve our modelling of many aspects 
of neutron star physics such as the effect of general rela- 
tivity on the r-modes, cooling rates, viscosity coefficients, 
magnetic fields, potenti al superfluidity, the formation of a 
solid crust etcetera (see Andersson & Kokkotas (2000 ) for 
a description of recent progress in these various directions) . 

In the following we will mainly consider uniform den- 
sity stars, i.e. use the gravitational-wave timescale given 
by (!]). In estimating the dissipation timescale t v u due 
to the presence of the crust, we need to evaluate t v u 
—2E/(dE/dt) where E is the mode-energy, and dE/dt fol- 
lows from an integral over the surface area at the crust- 
core boundary (assumed to be located at rad ius Rb), cf. 
equation (3) of Bildsten & Ushomirsky (2000). To evalu- 
ate this integral we use the standard result for the shear 
viscosity in a normal fluid. To incorporate our uniform 
density model, we make the reasonable assumption that 
the density of the star is constant (~ M/R 3 ) inside radius 
R. Then it falls off rapidly in such a way that the base 
of the crust is located at a radius only slightly larger than 
R. Hence, it makes sense to use Rb « R. If we neglect the 
small mass located outside radius R we can then immedi- 
ately compare the res ult for the viscous boundary layer to 
the tim escales used by Andersson, Kokkotas & Stergioulas 
(199S| ). In the end, our estimate for the dissipation due to 



the presence of the viscous boundary layer is 
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(2) 

which is a factor of 2 larger than that of Bildsten and 
Ushomirsky. This difference arises simply because the 
mode-energy E is this factor larger for uniform density 
models. The star is assumed to have a uniform tempera- 
ture distribution, with core temperature T . 

In Figure |l| we show the instability window obtained 
from our revised estimate. As is clear from this Figure, 
the presence of a viscous boundary layer in an old, rela- 
tively cold neutron star is, indeed, important. However, 
Bildsten and Ushomirsky's conclusion that the r-mode in- 
stability is irrelevant for the LMXB cannot be drawn from 
Figure [y. On the contrary, the Figure suggests that the 
instability may well be limiting the rotation of these sys- 
tems. 

3. THERMAL RUNAWAY IN RAPIDLY SPINNING NEUTRON 

STARS 

The fact that our revised instability curve for r-modes 
damped by dissipation in a viscous boundary layer agrees 
well with the fastest observed neutron star spin frequen- 
cies, cf. Figure EL motivates us to speculate further on the 
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relevance of the instability. We want to model how the 
potential presence of an unstable r-mode affects the spin- 
evolution of rapidly spinning, accreting neutron stars. To 
do this we use the phenomenol ogical two-parameter model 
devised by Owen ct al. (1998| ), which is centered on evo- 
lution equations for the rotation frequency Q and the (di- 
mensionless) r-mode amplitude a. Complete details of our 
particular version of this model will be given elsewhere. 

At the qualitative level, our results are not surprising. 
Accreting stars in the LMXB are expecte d to have core 
temp eratu res in the range 1 — 4 x 10 8 K (Brown & Bild- 
sten, ) . For such temperatures the dissipation due to 



the viscous boundary layer gets weaker as the tempera- 
ture increases. Consequently, the si tuation here is essen- 
tially identica l to t hat considered by Levin (1999|) (see also 
Spruit (199§1 ) and |Bildsten fc Ushomirsky (2000| )). After 



accreting and spinning up for something like 10 years, 
the star reaches the period at which the r-mode instabil- 
ity sets in. For our particular estimates this corresponds 
to a period of 1.5 ms (at a core temperature of 10 8 K). 
It is notable that this value is close to the 1.56 ms pe- 
riod of PSR1937+21. Once the r-mode becomes unstable 
(point A in Figure |l|) , viscous heating (now mainly due to 
the energy released in the viscous boundary layer) rapidly 
heats the star up to a few times 10 9 K. The r-mode am- 
plitude increases until it reaches a prescribed saturation 
level (amplitude a s ) at which unspecified nonlinear effects 
halt further growth (point B in Figure [j]) . Once the mode 
has saturated, the neutron star rapidly spins down as ex- 
cess angular momentum is radiated as gravitational waves. 
When the star has spun down to the point where the mode 
again becomes stable (point C in Figure |l|) , the amplitude 
starts to decay and the mode plays no further role in the 
spin evolution of the star (point D in Figure [l]) unless the 
star is again spun up to the instability limit. Two exam- 
ples of such r-mode cycles (corresponding to a s = 0.1 and 
1, respectively) are shown in Figure [fj. 

The real surprise here concerns the quantitative pre- 
dictions of our model. As already mentioned, the model 
suggests that an accreting star will not spin up beyond 
1.5 ms. This value obviously depends on the chosen stel- 
lar model, but it is independent of the r-mode saturation 
amplitude and only weakly dependent on the accretion 
rate (through a slight change in core temperature). In 
fact, the accretion rate only affects the time it takes the 
star to complete one full cycle. As soon as the mode be- 
comes unstable the spin-evolution is dominated by gravi- 
tational radiation and viscous heating. Once the star has 
gone through the brief phase when the r-mode is active it 
has spun down to a period in the range 2.8-4.8 ms (cor- 
responding to 0.01 < a s < 1). Based on these results we 
propose the following spin-evolution scenario: An accret- 
ing neutron star will never spin up beyond (say) 1.5 ms. 
Once it has reached this level the r-mode instability sets 
in and spins the star down to a several ms period. At this 
point the mode is again stable and continued accretion 
may resume to spin the star up. Since the star must ac- 
crete roughly 0.1M Q to reach the instability point, and the 
LMXB companions have masses in the range 0.1 — O.4M , 
it can pass through several "r-mode cycles" during its life- 
time. 

Let us confront this simple model with current obser- 
vations. To do this we note that our model leads to one 



main prediction: Once an accreting neutron star has been 
spun up beyond (say) 5 ms it must remain in the rather 
narrow range of periods 1.5 — 5 ms until it has stopped 
accreting and magnetic dipole braking eventually slows it 
down. Since a given star can go through several r-mode 
cycles before accretion is halted one would expect most 
neutron stars in LMXB and the MSP to be found in the 
predicted range of rotation rates. As is clear from Figure [l| 
this prediction agrees well with the range of rotation peri- 
ods inferred from observed kHz quasiperiodic oscillations 
in LMXB. The observed range shown in Figure [l] corre- 
spon ds to rotation frequ encies in the range 260-590 Hz 
(cf. [van dcr Klis (2000| )). Our model also agrees with 
the observed data for MSP, which arc mainly found in the 
range 1.56 — 6 ms, see Figure |[ In other words, our pro- 
posed model is in agreement with current observed data 
for rapidly rotating neutron stars. 

Finally, it is worthwhile discussing briefly the detectabil- 
ity of the gravitational waves that are radiated during the 
relatively short time when the r-mo de is saturate d and the 
star spins down. As was argued by Levin (1999), the fact 
that the r-mode is active only for a small fraction of the 
lifetime of the system (something like 1 month out of the 
10 7 years it takes to complete one full cycle) means that 
even though these sources would be supremely detectable 
from within our galaxy the event rate is far too low to make 
them relevant. However, it is interesting to note that the 
spin-evolution is rather similar to that of a hot young neu- 
tron star once the r-mode has reached its saturation am- 
plitude. This means that we can analyze the detectability 
of the emerging grav itational waves using the framework 
of pwen et al. (1998 ) . We then find that these events can 
be observed from rather distant galaxies. For a source in 
the Virgo cluster (assumed to be at a distance of 15 Mpc) 
these gravitational waves could be detected with a signal 
to noise ratio of a few using LIGO II. However, even at the 
distance of the Virgo cluster these events would be quite 
rare. By combining a birth rate for LMXB of 7 x 10 -6 
per year per galaxy with the fact that the the volume of 
space out to the Virgo cluster contains ~ 10 3 galaxies, and 
the possibility that each LMXB passes through (say) four 
r-mode cycles during its lifetime we deduce that one can 
only hope to see a few events per century in Virgo. In 
order to see several events per year the detector must be 
sensitive enough to detect these gravitational waves from 
(say) 150 Mpc. This would require a more advanced de- 
tector configuration such as a narrow-banded LIGO II. We 
will discuss this issue in more detail elsewhere. 

4. ADDITIONAL REMARKS 

Before concluding our discussion, we recall that the ini- 
tial excitement over the r-mode instability was related to 
the fact that it provided an explanation for the relatively 
slow inferred spin rates for young pulsars. In view of this, 
it is natural to digress somewhat and discuss how the pic- 
ture of the r-mode instability in hot, newly born neutron 
stars is affected by the possible formation of a solid crust. 
Hence, we consider the evolution of a neutron star just af- 
ter its birth in a supernova explosion. At a first glance, we 
might expect to model its r-mo de amplitude in the stan- 
dard way, cf. |Owenet al. (1998[ ), using the normal (crust- 
free) fluid viscous damping times for stellar temperatures 
above the melting temperature of the crust (T m ), and the 
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Fig. 1. — The r-mode instability window relevant for old neutron 
stars. We show results for the simplest (crust-free) model, where 
gravitational radiation reaction is balanced by standard shear vis- 
cosity at low temperatures and bulk viscosity at high temperatures 
(thin solid line). Also shown are Bildsten and Ushomirsky's estimate 
for a star with a crust (dashed line) and our improved estimate of 
this situation (thick solid line). The r-mode is potentially unstable 
in the region above each curve. The illustrated results correspond 
to a canonical neutron star for which mass shedding at the equator 
sets in at the Kepler period Pk ~ 0.8 ms. We illustrate two typi- 
cal r-mode cycles (for mode saturation amplitudes a a = 0.1 and 1), 
resulting from thermo-gravitational runaway after the onset of in- 
stability. Once accretion has spun the star up to the critical period 
(along the indicated spin-up line (thick vertical line)) the r-mode 
becomes unstable and the star evolves along the path A-B-C-D. Af- 
ter a month or so, the mode is stable and the star will cool down 
until it again reaches the spin-up line. For comparison with obser- 
vational data, we indicate the possible range of spin-periods inferred 
from current LMXB data (shaded box) as well as the nh«prwH peri 
ods an d estimated upper limi ts of the temperature (cf. [Andersson 

>idly 



Kokko as & Stergioulas (1999)) of some of the most rapidly spinning 
MSP (short arrows) 
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We compare the "r-mode cycle" predicted for ac- 
neutron stars to the observed rapidly spinning pulsars. 



ho ohscrvcd data is taken from the Princeton 



pulsar database 

([ittp:/ /pulsar. princcton.edu/pulsar/catalog sht.m ) as well as a sam- 



ple ol MSF recently discovered in 4YTuc (Uamilo et aL, 1999). 



visco us boundary layer damping time for temperatures be- 



low T 7 n- However, the situation is a little more complicated 



than this. Recall that the latent heat (i.e. the Coulomb 
binding energy) of a typical crust is Ei at ~ 10 ergs, while 
the r-mode energy is E m « 2a 2 (l ms/P) 2 x 10 51 ergs. 
Provided that the time taken for the star to cool to T m 
is sufficiently long, the energy in the r-mode (which grows 
exponentially on a timescale t gw ~ 20 s) will exceed Ei at , 
preventing the formation of the crust, even when T < T m . 
Th en the star will spi n down in the manner described by, 
eg. |Owen ct al. (1998 ). This phase will end either because 
the star leaves the instability region of the f2 — T plot (see, 
for example, Fig. 1 of Owen et al.), or because the mode 
energy in the outer layers of the star (where the crust is 
going to form) has fallen below the crustal binding energy. 
We can estimate that this would happen at a frequency 
~ 70 Hz/a s , by equating Ei at to roughly 10 % of E m ). 
A more accurate treatment would take into account the 
local kinetic energy of fluid elements and since the latter 
is smaller near the poles than near the equator, the crust 
might form earlier at the poles. Clearly the problem of 
crust formation in an oscillating star requires further in- 
vestigation. The final spin period will be around 15 ms, 
if the r-mode grows to an amplitude of a s ~ 1, consistent 
with the extrapolated initial spin rates of many young pul- 
sars. On the other hand, if the mode is not given time to 
grow very large it will not prevent crust formation at T m . 



Such a scenario was described by Bildsten & Ushomirskv 
(200C| ), who noted that the r-mode instability would then 



not spin the star down beyond a much higher frequency 
Using our estimated timescales the resultant spin period 
would be 3 — 5 ms in this scenario. Which scenario ap- 
plies depends sensitively on the early cooling of the star, 
the crustal formation temperature and perhaps most im- 
portantly the initial amplitude of the r-mode following the 
collapse. It is, in fact, possible that both routes are viable 
and that a bimodal distribution of initial spin periods re- 
sults. A likely key parameter is whether the supernova 
collapse leads to a large initial r-mode amplitude ~ a s 
or not. An initial period of ~ 15 ms would fit the long 
established data for the Crab, whil e the recently discoy- 
ered 16 ms pulsar PRS J0537-6910 (|Marshall et al., 1998|) 
requires a considerably shorter initial period of a few ms. 

In conclusion, we have reexamined the effect that the 
dissipation due to a possible viscous boundary layer in a 
neutron star with a solid crust has on the stability of the r- 
modes. By combining our new estim ates with the the rmal 
runaway introduced by Levin (199S ) and Spruit (1998 ), we 
arrive at a spin-evolution model that agrees with present 
observations for rapidly spinning neutron stars. In partic- 
ular, our predictions agree well with observations of both 
LMXB and MSP. Furthermore, the model can potentially 
explain the extrapolated spin periods of the young pul- 
sars. Since it brings out this unified picture, our simple 
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model has many attractive features, and we are currently 
investigating it in more detail. 
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